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Abstract 

We investigate the issues of future oscillations around the phantom divide for f(R) gravity. 
For this purpose, we introduce two types of energy density and pressure arisen from the f(R)- 
higher order curvature terms. One has the conventional energy density and pressure even in the 
beginning of the Jordan frame, whose continuity equation provides the native equation of state 
Wde- On the other hand, the other has the different forms of energy density and pressure which 
do not obviously satisfy the continuity equation. This needs to introduce the effective equation of 
state w e ff to describe the /(i?)-fluid, in addition to the native equation of state w)de- We confirm 
that future oscillations around the phantom divide occur in f(R) gravities by introducing two types 
of equations of state. Finally, we point out that the singularity appears ar x = x c because the 
stability condition of f{R) gravity violates. 
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I. INTRODUCTION 



Supernova type Ia(SUN la) observations has shown that our universe is accelerating 



Also cosmic microwave background radiation [2j, large scale structure 



ing {4 1 have indicated that the universe has been undergoing an accelerating phase since 
the recent past. The standard model of ACDM enables to explain these observational re- 
sults within observational error bound. However, this model suffers from the cosmological 
constant problem and thus, one needs to find another model. Up to now, the /(-R)-gravity 
as a modified gravity model remains a promising model to explain the present accelerating 
universe {5- 9]. f(R) gravities can be considered as Einstein gravity (massless graviton) with 
an additional scalar. For example, it was shown that the metric-/ (R) gravity is equivalent 
to the cjbd = Brans-Dicke (BD) theory with the potential [9|. Although the equivalence 
principle test in the solar system imposes a strong constraint on f(R) gravities, they may 
not be automatically ruled out if the Chameleon mechanism is introduced to resolve it in the 
Einstein frame. It was shown that the equivalence principle test allows f(R) gravity mode 



3j, and weak lens- 



that are indistinguishable from the ACDM model in the background universe evolution [10]. 

In order to point out the difference between ACDM and f(R) gravity, it is necessary to 
introduce the equation of state parameter u>de- Working with the /(-R)-gravity action in 



the Jordan frame 



11| . one has to use the different energy density and pressure (pde,Pde) 



n 



in compared to (pde;Pde) i n the Einstein-like frame Jl2l . Il3j |. This corresponds to the 
scalar-tensor (Brans-Dicke) theory in the Jordan frame 14j. In the Einstein-like frame, one 
needs only the native equation of state wde = Pde/pde as in the scalar-tensor (quintessence 
model) theory, while one requires two equations of states: Wbe = Pde/pde and the effective 



equation of state w e s 



151 ] because of non-minimal coupling of scalar to the gravity. It is 



worth noting that there is an essential difference between Einstein-like and Einstein frames 
because the latter is recovered from the conformal transformation in Jordan frame 



arified that although there is 

a phantom phase when using the native equation of state 17(, there is no phantom phase 



For the holographic dark energy model, two of authors have c 



when using the effective equation of state [l8|. 
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Recently, there were a few of important works which explain the oscillation around the 
future de Sitter solution with w^s = — 1 using /(i?)-gravity [3] and its Brans-Dicke the- 



ory 



141 ] . Interestingly, the authors in 13| have shown that the number of phantom divide 



2 



crossings are infinite when using the Ricci scalar perturbation, which is confirmed by ana- 
lytical condition and numerical way. 

In this work, we focus on the issues of future oscillations around the phantom divide for 
f(R) gravity. In order to confirm the appearance of future oscillations around the phantom 
divide w^s = — 1, we introduce two types of equation of states wde(^de) and w e s arisen 
from the /(-R)-fluid. We clarify the difference between two different sets of energy density 
and pressure by observing the "negative and effective" equations of state. Finally, we point 
out that the singularity appears at x = x c because the stability condition of f(R) gravity 
violates when F' — f" — at the certain point x = x c . 



II. FUTURE VOLUTION WITH /(fl)-FLUID IN EINSTEIN-LIKE FRAME 

We start from the action of f(R) gravity with matter as 

1 = 2^ / ^^fW + UWm), (1) 

where f(R) is a function of Ricci scalar R with n 2 = 8irG and 7 m is the action for matter 
which is assumed to be minimally coupled to gravity. Here the action / is initially written 
in Jordan frame and i[) m denotes matters. Taking the variation of the action ([T]) with respect 
to metric g^ u , one obtains 

FGfiv = k 2 T^ - l -g, u (FR - f) + V M V V F - g, u V 2 F, (2) 

where G^ u = R^ v — ^-R is the Einstein tensor and F(R) = f'(R). Assuming the flat 
Friedmann-Roberston- Walker (FRW) universe 



dsl RW = -dt 2 + a 2 it) [dr 2 + r 2 dttf) (3) 



with a(t) is the scale factor, we obtain the two Friedmann equations from (J2]): 

3FH 2 = k 2 P m + l -{FR-f)- 3HF, (4) 

-2FH = k 2 (p M + Pm) + F - HF, (5) 

where H = a/a is the Hubble parameter, the overdot denotes the derivative with respect to 
the cosmic time t, and pm and pm are the energy density and pressure of all perfect fluid-type 
matter, respectively. On the other hand, the scalar curvature R defined by 

R = q(h + 2H 2 ) (6) 



plays an independent role in the cosmological evolution because we are working with f(R)- 
fluid. For our purpose, we introduce the new variable x = In a, then (j3J) and ([6]) take the 
forms 



H 



R 



6 H 



dH 

dx 



dx 



2H' 



_(f-m-H 2 F'—-< Pm 

6 dx 



K 2 p® m e 4x 



X, (7) 
(8) 



where p® m is the current density of cold dark matter (CDM) and % = Pr/Pm — 3.1 x 10 -4 
is the current density ratio of radiation and dark matter. Regarding (J7|) as the evolution 
equation, we rewrite it clS db compact form 



K 



H 2 = y (pde + Pn 



(9) 



where pde, Pm, and p r represent dark energy, dark matter and radiation, respectively. 
Comparing ([7| with (Q leads to a definition of dark energy density arisen from the f{R)- 



gravity 



12 



13j 



Pde — — 



(F — 1) I H^- + H 2 )--(f-R)- H 2 F'^ 



dx 



6 



dx 



This dark energy density satisfies the conservation law as 

Pde + 3H(l + w D e) Pde = 
with the native equation of state 



Pde 
Pde 



(10) 



(11) 



(12) 



Importantly, we observe that even starting with the f(R) gravity action in the Jordan frame, 
we have manipulated it so that the Einstein equation <^ is rewritten effectively as 



G 



G^(l — F) + k 2 T^ - 1^ (FR -f) + V M V,F - g, v V 2 F, 



(13) 



to derive <^ and for obtaining the standard dark energy pde and pressure pde- Hence we 
call the solution to OH]) with ffTUj) as cosmological evolution with f(R) -fluid in Einstein-like 
frame. 

In order to solve Q and (jSj) simultaneously, we define the convenient variables call the 
reduced Ricci scalar r and the present matter-density parameter 



r = 



Pcrit 



3H 2 ' 



(14) 



and density parameters 

fi m = — , f2 r = , ^DE = 1 — fim — fir- (15) 

Pcrit Pcrit 

Then two equations (JZj) and (jSJ) can be written four equations in terms of new variables 

-(2r-l)^ m , (16) 



dx 

§ = (17) 

f = (r-2)tf. (19) 



Considering 



^ = _ 2 „ DE(r _ 2) + _y^, (20) 

one obtains the native equation of states as functions of r and density parameters as 

, n s . 2r-4 + 3fi m + 4fi r 
WDE (r, n m , fi r ) = -1 - 3(1 _ nm _ a) • (21) 

At this stage, we wish to comment on the tilde-definition for density parameters used in 
Ref. 



12J 



fim — o — fi m /i , (22) 

Pcrit 

f2 r = f^L = Q v h 2 , (23) 

Pcrit 

fi DE = 4^ = fi D E /i 2 , (24) 

Pcrit 

which are clearly different from our definition by h 2 factor, which is defined as h(x) = 
H(x)/Hq. Finally, we mention that the initial conditions for Q m ,Q v ,r and H are given by 

fi m (0) = <, a(0) = xfim, r(0) = 2 + /n, if(0) = l, (25) 

where /ii is the time derivative of H at the present time t = with a(0) = 1. Then, 
x G [—00, 00] with x = at the present time. We note that hi is related to deceleration 
parameter q as 

h 1 = -(l + q) (26) 
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with 



aa 

q = - — 

GT 



a 1 



1 dH 



H dx 



1 



-fl! ~ 1. 



(27) 



Now we study the cosmological evolution by choosing four specific models of /(-R)-gravity 
with the native equation of state wde only. 



A. Cosmological constant 

For cosmological constant case, the function f(R) is simply given by 

f(R) = R + A. (28) 

In this case one cannot use ([IS]) because of F'(R) = f"(R) = 0, but other equations are 
being used to derive the solutions. Equation (JTj) becomes 

H 2 = -± + H 2 n° m (e- 3 * + X e- 4x ). (29) 
Differentiating this equation with respect to x, we get 

2H^- = H*i& (-3e~ 3 * - Axe-**) . (30) 
Plugging this into the definition of r, one gets 

r = 2 - - 2Q r . (31) 
Hence, the relevant equation is just 

(IT i 

— = -(2r-l)fi m + 4ra (32) 
ax 2 

with j«de = — 1. Fig. [I] depicts the evolution for the cosmological constant like the ACDM 
without future oscillations around the phantom divide %e = ^as — — 1 • We point out 
that the reduced Ricci scalar r takes the value of r^s = 2 because of Rds = 12H 2 for the de 
Sitter spacetimes. 

B. Power-law gravity 

When f(R) takes the power- law form 

f(R) = R + f R a , (33) 
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% 




X 



FIG. 1: Time evolution of cosmological parameters for cosmological constant case: ^DE(blue), 
^DE(green), f2 m (magenta), O r (brown), and r (yellow) for fi^ = 0.23 and \ — 3.04 x 10~ 4 ,/ii = 

Its derivatives with respect to R are given by 

F(R) = 1 + af R a ~\ F'(R) = a (a - l)f R a ~ 2 - (34) 

wde and r in Figs. |2] and [3] represent the future oscillations around the phantom divide for 
a = 1/2 and 1/3, respectively but they do not show past oscillations around the phantom 
divide. 
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FIG. 2: Time evolution of cosmological parameters for power-law case: ^DE(blue), f^DEfereen), 
O m (magenta) , f2 r (brown), r (yellow), and F' (black) for /o = —3.6, a = 1/2,0^ = 0.23, \ = 
3.04 x 10- 4 ,hi = -(I + f)Q^. 

C. Exponential gravity 

Now we wish to apply the result of previous sub-sections to an exponential gravity. Firstly, 
when the function f(R) is given by 

f(R) = R-(3R s (l-e- R/R °), (35) 
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FIG. 3: Time evolution of cosmological parameters for power-law case: tt>DE(blue), firjE (green) , 
O m (magenta), f2 r (brown), r (yellow), and F' (black) for /o = —6.5, a = 1/3,0^ = 0.23, \ = 
3.04 x 10" 4 ,/ii = -(I + X)Sl° m . 



its derivatives with respect to R are given by 

-R/R s p>s R \ _ P_ p -R/R s 



F(R) = 1-Pe 



(36) 



Fig. H] indicates no the appearance of future (past) oscillations around the phantom divide 
for a given parameter R s = —0.05 and (3 = 1.1 which is the same results found in 20(. We 
note that F' does not appear in Fig. 4 because its value is extremely large as 10 87 . 
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FIG. 4: Time evolution of cosmological parameters for an exponential case: u>DE(blue), Ode (green), 
f2 m (magenta), f2 r (brown), r (yellow), and F' (black) for R s = — 0.05, /3 = 1.1, = 0.23, \ = 
3.04 x 10" 4 ,/ii = -(I + f)Q^. 

D. Hu and Sawicki model 



The Hu and Sawicki model takes the form 

R 2 ' 



f(R) = R-/2R c 



1 



Rl 



(37) 
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Its derivatives with respect to R are given by 

-(n+l) 

F(R) 



R ( R 2 



(38) 



F\R) 



2 fin 

Rr 



1 - (2n + 1 



R 2 



1 



(39) 



R 2\ -(n+2) 

Fig. [5] shows that future oscillations around the phantom divide was = — 1 appears for 
i? c = —1.0, fi = —1.5, and n = 2, but there is no evolution toward the past direction from 
x c = 0. Similarly, there are future oscillations around r^s = 2 for the reduced Ricci scalar r. 
We will explain why the evolution to the past is not allowed in the Hu and Sawicki model 
in section IV. 



III. COSMOLOGICAL EVOLUTION IN JORDAN FRAME 



In this section we derive evolution equations in Jordan frame without manipulation. From 
equations (jl]) and (J5J), we have 



H 2 



-H 



3F 
1 

2F 



K 2 p M + -(FR- f)-3HF 



K 2 {p M +p M ) + F-HF 



Introducing the reduced Ricci scalar r 
R 



QH 2 ' 



we can obtain an important relation 

F ' dR i 1 

F^ = - 1 + Q - + n < + 6H 2 



D / 

~F 



(40) 
(41) 

(42) 
(43) 



In this case, we read off energy density and pressure of the dark energy component from 
gO]) and flUD as [11] 

1 (FR-f 



Pde 



Pde 



k 2 F I 2 

FR-f 



k 2 F { 2 
Rewriting the Einstein equation (j2J) as 

^(m) 
2 1 V U , ,.2/tiDE 



3HF 

f 2HF + Fy 



+ AC T. 



(44) 
(45) 

(46) 
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FIG. 5: Time evolution of cosmological parameters for the Hu and Sawicki model: 'WDE(blue), 
^DEfereen), fi m (magenta), J7 r (brown), r (yellow), and F' (black) for R c = —1.0,fi = — 1.5, n = 
2, n° m = 0.23, X = 3.04 x 10~ 4 , hi = -(§ + § )n° m . 



we obtain the non-conservation of continuity relation by requiring the Bianchi idensity 



Pde + (p DE + p DE ) = j^Pm- 



(47) 



Because of the non-zero coupling term between pu and J^, we must define an "effective" 
equation of state for /(i?)-fluid 

-1 + O m + O r + r — Qjfip 



(48) 
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with the native equation of state u>de = Pde/pbe- This is similar to the Brans-Dicke theory 
approach [lij] . Defining density parameters newly as 



relevant quantities are expressed by 

R = QH\ (50) 
3H 2 

PDE = —5- {1 - - ^r} , (51) 



Pde = {-1 - - |(r-2) [> . (52) 



Insereting these into the definition of native equation of state, one gets 

WDE = - 1 " 3(1 -n m - a) ' (53) 

which is the same as Eq. (1211) but their forms of f2 m and f2 m are different from those in (|2ip . 
Four equations to be solved become 

^2r-2 + Q m + n r + r- -J—] fi m , (54) 



dx V 6# 2 F 

^2r - 1 + fi m + a + r - — ) n T , 



dx V 6# 2 F 

dr F ( f \ 

S = - 2r(r - 2) + 67^ (- 1 + °- + + 6/M • < 56) 

Finally, we have to consider the initial conditions 

H(0) =tfo = l, a(0) = X^m, r = 2 + /n = -(1 + q), (58) 

with g the deceleration parameter. 

Now we solve the above four differential equations together with initial conditions by 
selecting four interesting models. 

A. f(R) = R + A case 

In this case, we cannot use Eq. ( 156]) . which is valid only for F' ^ 0. In this case, r is 
given by 

r = 2 - -n m - 2Q T . (59) 
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Hence the density and pressure of /(-R)-fluid are given 
3H 2 

pde = —2- (i - a* - a) , (60) 

3H 2 

pde = — (-i + a n + a) • (6i) 

Therefore, we have wde = — 1 = w e s- 

^ = -(2r-l)fi mi (62) 

^a ~ 

— - = -2ra, 63 
ax 

r = --fi m - 2a, (64) 

§ = (r-2)2T. (65) 

Note that these equations are exactly the same as the previous section. Fig. [6] shows the 
same result that the future oscillations around the phantom divide does not appear in the 
Jordan frame. We observe that the deceleration parameter q is fixed by the relation 

1 + <? = -Q + £K- (66) 

For n DE = 0.77, n° m = 0.23, x = 3.1 x 10~ 4 , it gives us q = -0.6549643500 and h x = 
-0.3451398400. 

B. f{R) = R + f R a case 

In this case, F(R) is given 

F(R) = 1 + af R a ~\ (67) 

Its derivatives with respect to R are given by 

F(R) = l + af R a ~\ (68) 
F'(R) = a(a - l)hR a ~ 2 - (69) 

Figs. [7] and [8] show future oscillations around the phantom divide using wde as in Figs. 2 
and 3 but the past evolution is terminated near x c ~ —2.4 for a = 1/2 and x c ~ 3.4 for 
a = 1/3. Also, the reduced Ricci scalar r shows similar oscillating behaviors. This indicates 
a violation of stability condition as will explain in section IV. We confirm the appearance 
of future oscillations around the phantom using the effective equation of state w e Q. 
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FIG. 6: Time evolution of cosmological parameters for cosmological constant case, iude (blue) , 
^DE(green), Q m (magenta), f2 r (brown), r (yellow), for fi^ = 0.23 and x = 3-04 x 10~ 4 , hi = — (§ + 

C. Exponential gravity 

For the exponential gravity, the function f(R) is given 



f(R) = R-(3R s (l- e -W). 



(70) 
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FIG. 7: Time evolution of cosmological parameters for power-law case in Jordan frame: ?Dde (blue) , 
io e g(cyan), ^DEfere-en), f2 m (magenta), f2 r (brown), r (yellow), and F' (black) for /o = —3.6, a = 
1/2,^ = 0.23, X = 3.04 x lO" 4 ,^ = -(§ + %)tt° m . 



Its derivatives with respect to R are given by 
F(R) = l-f3e~ RIR % 



F\R) = 

K s 



-R/R s 



(71) 
(72) 



Figs. |9] depicts that future oscillations around the phantom divide does not appear when 
using wbe and w e g. There is no essential difference between Einstein-like (Fig. 4) and 
Jordan frames (Fig. 9). We would like to mention that F' does not appear in Fig. 4 because 
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FIG. 8: Time evolution of cosmological parameters for power-law case in Jordan frame: ?Dde (blue) , 
w e R (cyan) , ^DEfereen), f2 m (magenta), f2 r (brown), r(yellow), and F' (black) for /o = —6.5, a = 
1/3,0°, = 0.23, x = 3.04 x lO" 4 ,^ = -(§ + f)^. 

its value is extremely large as 10 87 . 
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FIG. 9: Time evolution of cosmological parameters for exponential case in Jordan frame: 
?Z>de (blue) , w e g (cyan) , Ode (green), (magenta), O r (brown), r (yellow), and F' (black) for 
R s = -0.05, [3 = 1.1, n% = 0.23, x = 3.04 x 10" 4 , hi = -(§ + §)^. 



D. Hu and Sawicki case f(R) = R — fiR c 



i-u+S 



Its derivatives with respect to R are given by 

R ( R 2 ^ ~^ n+1 ' ) 



F'(R) 



2 fin 

Rr 



1 - (2n + 1 



Rl 



R^_ 
Rl 



-(n+2) 



(73) 
(74) 
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FIG. 10: Time evolution of cosmological parameters for Hu and Sawicki case in Jordan frame: 
WDE(blue), w e g (cyan) , f^DE (green), Cl m (magenta), Q T (brown), r (yellow), and F' (black) for R c = 
-1.0, /i = -1.5, n = 2, = 0.23, x = 3.04 x 10" 4 , hi = -(§ + | )0^. 

Fig. [TUI shows that future oscillations around the phantom divide w^s = — 1 appears for 
R c = —1.0, = —1.5, and n = 2 using two equations of state wbe and ty e fr but there is no 
evolution toward the past direction from x c = 0. Also, there are future oscillations around 
rds = 2 for the reduced Ricci scalar r. This will be explained in the next section. This 
confirms the results (Fig. 5) in the Einstein-like frame. 
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IV. SINGULARITY IN COSMOLOGICAL EVOLUTIONS 



In this section, we wish to explain the singularities encountered in the cosmological evolu- 



bounds: 



13| 



tion of /(i?)-fluid 21] . First of all, we mention that /(-R)-gravity should satisfy the following 



f'(R) = F(R) > 0, f"(R) = F'(R) > 0. (75) 

These are necessary to guarantee that the Newtonian gravity solutions are stable and that 
the matter-dominated stage remains an attractor with respect to an open set of neighboring 
cosmological solutions in f(R) -gravity. In the perturbation theory, the former is necessary 
to show that the gravity is attractive and and the graviton is not a ghost, whereas the latter 
needs to ensure that the scalaron of a massive curvature scalar does not have a tachyon. 
We show how the singularities appear from the /(i?)-fluid. From the observation of two 
equations (|18p and (]56|) which are equivalent to two first Friedmann equations, the second 
term involves F'(R) as the denominator. Hence, if F'(R) = at a certain point of x = x c , 
it gives rise to an singularity at which some cosmological parameters blow up. This results 
from the violation of stability condition f(R) gravity: the second of f!75|) . 

In order to show the presence of singularities explicitly, we use the graph of F'( 
function of x. From Fig. 5, we find that the singularity appears at x c = and thus, the 
backward evolution is not allowed in Einstein-like frame. In the Jordan frame, the power-law 
gravity shows the singularities at x c ~ —2.4 for a = 1/2 (see Fig. 7) and x c ~ —3.4 for 
a = 1/3 (see Fig. 8), while from Fig. 10, we find that the singularity appears at x c = and 
thus, the backward evolution is not allowed in Jordan frame. This completes the presence 
of singularities in the cosmological evolution of /(-R)-fluid. 



V. DISCUSSIONS 



We have investigated the issues of future oscillations around the phantom divide for f(R) 
gravity by introducing two types of energy density and pressure arisen from the /(-R)-fluid. 
One has the conventional energy density pde and pressure pde even in the beginning of 
the Jordan frame, whose continuity equation provides the native equation of state u>de = 
Pde/pde- Hence, we call this frame as the Einstein-like frame. 
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On the other hand, the other has the different forms of energy density pde and pressure 
Pde which do not obviously satisfy the continuity equation. This needs to introduce the 
effective equation of state w c q to describe the /(-R)-fluid precisely, in addition to the native 
equation of state wde = Pde/pde- We confirm that future oscillations around the phantom 
divide always occur in /(i?) -gravities by introducing two types of /(-R)-gravity: one is the 
power-law potential (133]) with the exponent a = 1/2 and 1/3 and the other is the Hu and 
Sawiciki model fl37|) . In the Jordan frame, the former did not show past oscillations around 
the phantom divide and its evolution was terminated around x c ~ —2.4 for a = 1/2. On 
the the other hand, the latter did not provide any past evolution in both Einstein-like and 
Jordan frames. Similarly, we confirm that there are future oscillations around r^s = 2 for 



the reduced Ricci scalar r 



13]. As was expected, the cosmological constant model has no 



frame- dependence and we could not find any future oscillations around the phantom divide 
around u>de = — 1 for the exponential gravity in (1351) . 

For whole evolution from the past to future when imposing initial conditions at the 
present time, the cosmological evolution is allowed in the Einstein-like frame better than in 
the Jordan frame. This means that the cosmological evolution of /(-R)-fluid determined from 
its form of energy density and pressure, depending on the given frame. Also, it was proven 
that the termination (singularity) appeared in cosmological evolution is closely related to 
the form of /(-R)-fluid for given frame. This has arisen from F' = in the first Friedmann 
equations (Tl8]) and ( )56l) . As a result, it is so because of the violation of the stability condition 
(non-tachyon) of f{R) gravity. 

Consequently, we have successfully performed (whole) cosmological evolution of f(R) 
gravities by choosing two different state variables of energy density and pressure, and pointed 
out why the singularity appeared in the backward evolution when the initial condition was 
chosen as the present time. 
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